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In this paper we extend both standard fault tolerance theory and Kitaev's model for quantum 
computation, combining them so as to yield quantitative results that reveal the interplay between 
the two. Our analysis establishes a methodology that allows us to quantitatively determine design 
parameters for a quantum computer, the values of which ensure that an overall computation of in- 
terest yields a correct final result with some prescribed probability of success, as opposed to merely 
ensuring that the desired final quantum state is obtained. As a specific example of the practical 
application of our approach, we explicitly calculate the number of levels of error correction con- 
catenation needed to achieve a correct final result for the overall computation with some prescribed 
success probability. Since our methodology allows one to determine parameters required in order to 
achieve the correct final result for the overall quantum computation, as opposed to merely ensuring 
that the desired final quantum state is produced, our method enables the determination of com- 
plete quantum computational resource requirements associated to the actual solution of practical 
problems. 

PACS numbers: 03.67.Lx, 03.67.Pp 
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I. INTRODUCTION 

The purpose of any computation, whether imple- 
mented by a quantum computer or a classical computer, 
is to compute the value of a function for specified values 
of the variables on which the function depends. In the 
case of a quantum computer, the final result is obtained 
from the outcome of a final measurement performed on 
the final quantum state produced by the quantum com- 
putation. This is illustrated schematically in Figure [TJ 
which also serves to fix the terminology we shall use in 
this paper to describe the quantum computation and the 
final measurement. In particular, in this paper, the term 
"quantum computation" very specifically refers solely to 
the dynamical evolution of the qubits from some initial 
quantum state to some final quantum state. Our use of 
this term does not include the subsequent, final measure- 
ment from which the final result of the overall computa- 
tion is obtained. We thus distinguish between the final 
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FIG. 1: Basic definitions. 
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quantum state of the quantum computation and the final 
result of the overall computation, of which the quantum 
computation is merely a part. 

Since the final measurement produces, in general, a 
probabilistically distributed set of outcomes, the ques- 
tion arises as to whether or not reliable final results for 
the overall computation can be obtained. An affirmative 
answer to this question is necessary in order that algo- 
rithms such as Shor's algorithm [1] can be successfully 
applied to practical problems. A partial answer to this 
question follows from the quantum computational model 
formalized by Kitaev ([2|, §4.1). Kitaev's model identifies 
a bound, p, on the probability that the final result of the 
overall computation is incorrect due to the indeterminism 
of the final measurement. It follows that reliable over- 
all computation can be achieved for sufficiently small p. 
However, Kitaev's model assumes that the quantum com- 
putation that produces the final quantum state (on which 
the final measurement is performed in order to give the 
final result of the overall computation) is perfectly imple- 
mented, with no errors. In other words, those problems 
that are within the purview of fault tolerance analysis, 
namely, the occurrence of errors in the quantum computa- 
tion, and the effectiveness of error correction techniques 
in reducing the effect of such errors on the final quantum 
state, are not addressed in Kitaev's model. There are 
thus two sources of error that may affect the final result 
of the overall computation: (1) errors which arise in the 
course of the quantum computation, and (2) errors due 
to the indeterminacy intrinsic to the final measurement 
which follows the quantum computation [3( . Fault tol- 
erance theory 0, @| addresses the first type of error, 
but does not consider the second. Kitaev's model ad- 
dresses the second type of error, but not the first. In the 
present paper, we extend both fault tolerance theory and 
Kitaev's model so as to take into account the combined 
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effects of both sources of error, including the interplay 
between them. 



In order to quantify the connection between fault tol- 
erance theory and Kitaev's model, we utilize a measure 
of the difference between the desired final quantum state 
(that would arise in the case of perfect quantum compu- 
tation) and the actual final quantum state (that is pro- 
duced in an actual quantum computation realized by a 
practical device). We define this measure in terms of 
a suitable norm of the difference between the two afore- 
mentioned states. We refer to this measure as the "imple- 
mentation inaccuracy." We then introduce a fundamental 
inequality, to which we refer as the Quantum Computer 
Condition (QCC), which requires that the implementa- 
tion inaccuracy be less than some prescribed bound. This 
fundamental inequality furnishes a criterion for successful 
fault tolerant quantum computation. We shall see that 
the use of this criterion allows the conditions that ensure 
fault tolerance to be incorporated into Kitaev's model in 
a remarkably straightforward way. 



Our results provide a quantitative relationship between 
fault tolerance theoretic constraints on the one hand, and 
the accuracy of the final result of the overall computation 
on the other. Consequently, our results can be applied to 
explicitly and directly determine requirements on practi- 
cal fault tolerant design parameters in order to achieve a 
specified accuracy for the final result of the overall com- 
putation. This goes beyond standard fault tolerance the- 
ory, which only addresses the accuracy with which the 
quantum computation realizes the desired final quantum 
state, and does not directly determine the accuracy of the 
final result of the overall computation. As an example of 
the application of our result, we show how to explicitly 
calculate the number of levels of concatenation of error 
correction required to directly achieve a specified proba- 
bility that the overall computation produces the correct 
final result, and not merely the desired final quantum 
state. 



II. THE KITAEV MODEL 

We begin by reviewing the Kitaev model @ , which can 
be depicted by the following diagram: 
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This diagram expresses the fact that the output of the 
quantum computation U is intended to be used in com- 
puting the function F, where F : X — » Y is an instance of 
the overall computational problem. Here I pU re is an ini- 
tialization map, which maps the input space X into pure 
states of the logical Hilbert space H logical) and O pur e is 
the corresponding readout map, which maps the output 
of the operation U onto the output space Y . In gen- 
eral Op Ur e is a measurement given by a projection- valued 
measure (or more generally a POVM) {E y } y£ y The 
symbol in the center of the diagram denotes the prob- 
abilistic inaccuracy, p, associated to the output of the 
overall computation. If p = 0, the computation always 
produces the correct result, and the diagram commutes. 
Thus, the quantity 1 — p furnishes a lower bound on the 
probability of the success of the overall computation. We 
note that p may be non-vanishing, even if the quantum 
algorithm is perfectly implemented, due to the fact that 
the measurement of the final quantum state produced by 
the quantum computation is necessarily quantum prob- 
abilistically distributed. Given an input x G X, the final 
result of the overall computation is distributed according 
the probability measure on Y as follows: the probability 
Pr x (y) of a singleton y G Y is (I pure (a:)|U' t E a U'|Ip Ur e(:E)}. 
Then, the near commutativity of Diagram Jl} means that 
for each x G X, the probability measure Pr x (y) is con- 
centrated at y = F(x), so that F{x) is the final result 
of the overall computation with high probability. Thus, 
Kitaev's formulation Q requires that the Diagram ([1]) be 
nearly commutative in a probabilistic sense, that is 



Or 



3 (a;)|[/ t E F(:E) [/|Ip Ure (x)) > 1 -p 



(2) 



This paper is organized as follows. In Section 2 we re- 
view the Kitaev model and re-express it in a form that 
is convenient for our analysis. In Section 3 we introduce 
the implementation inaccuracy and rigorously define the 
fundamental inequality known as the QCC. In Section 4 
we derive the constraints implied by combining the QCC 
with the Kitaev model. As an example of the practical 
utility of our general result, we explicitly calculate the 
number of levels of error correction concatenation needed 
to achieve a correct final result for the overall computa- 
tion with some prescribed success probability. In Section 
5 we present our conclusions. 



In other words, the measurement of the final quantum 
state gives the correct final result for the overall compu- 
tation with probability greater than 1 — p. If p is suf- 
ficiently small, e.g., p < 1/2 in the case of Y binary, 
a majority vote algorithm will successfully identify the 
correct outcome y = F(x). 

We wish to extend the Kitaev model to apply to the 
scenario in which errors may occur in the implementation 
of the quantum computation. As a preliminary step, we 
re-express the Kitaev model in a more convenient form 
for our purposes by generalizing Diagram §Q and eq.© 
through the replacement of wavefunctions with density 
matrices. This allows us to discuss the effect of errors 
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that transform pure states into mixed states. The gener- 
alized diagram is 



T( -ff logical) 
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where the action of the map G is given by G : p i— > UpW , 
p is the density matrix representing the initial state of 
the quantum computer, the action of the map Idm is 



given by l dr 



|Ipure(^))(Ipure(a;)|, O dm is the 



measurement corresponding to O pure , except that Odr 
acts on density matrices, and T(H) is the Banach space 
of trace class operators on a given Hilbert space, H. 
The probability Pr x (y) of a singleton y G Y is now 
tr(,yi^UIdm{x)U^ y/F^). The near commutativity of 



Diagram ([3]) means that 

tr(O dm (G (Idm (a;)))) 
tr 



^E F{x) Ul dm (x)U^E F{x) ) > 1-p (4) 



for all a; G X. Eq.Q generalizes eq.([2]) to states de- 
scribed by density matrices. 



III. INCLUSION OF RESIDUAL ERRORS AND 
THE QCC 

We now extend Kitaev's model by allowing for the in- 
evitable survival of residual errors in any realistic im- 
plementation of a quantum computation, even upon 
successful application of fault-tolerance techniques. In 
other words, we will extend the analysis to include fault- 
tolerant operation, but in such a way as to ensure a pre- 
scribed success probability to achieve the correct final 
result of the overall computation. 

We consider a function F, and an implementation of 
an overall computation (including a quantum computa- 
tion followed by a final measurement) that is intended 
to calculate F for some value of the input. The rela- 
tionship between the ideal, error-free quantum computa- 
tion, G(p) = U pU\ defined on the logical Hilbert space, 
-ffiogicai, and the actual dynamical map, P, implemented 



by the physical quantum computer on the computational 
Hilbert space, ff comp , can be described by the following 
diagram: 



T(-£/c 0m p) 
M { i_ c} 

T(i? logical) 
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T(_/? CO mp) 

M {c ^ 1} 

T(.ff logical) 



(5) 



where the superoperators M|i^ c } and M{ c ^i} are link- 
ing maps that are mathematically required to relate 
states in the logical space Pi og icai to states in the com- 
putational space H comp . This is because, in general, 

Plogical 7^ ^comp [jj] • 

We formally express the content of Diagram ([5]) via a 
relation between an idealized (G) and actual (P) quan- 
tum computation: 



\\M {c ^ l} (P ■ (M {l ^ c} (p))) - Giplh <a 



(0) 



where || • ||i signifies the Schatten 1-norm. The quantity 
on the left-hand-side of this inequality furnishes a mea- 
sure of the inaccuracy of the actual implementation of G 
by P. It tells us how well a practical quantum computing 
device actually implements an ideally defined quantum 
computation. We will refer to the left-hand-side of this 
inequality as the implementation inaccuracy of the quan- 
tum computation. The quantity, a, on the right-hand- 
side of this inequality is a bound on the implementation 
inaccuracy. The value of a is prescribed as a requirement 
on the performance of the quantum computation. Thus, 
the diagram states that P can implement the idealized, 
perfect quantum computation, G, with an inaccuracy no 
greater than a. We refer to the entire inequality as the 
Quantum Computer Condition (QCC). 

Thus, we have that Diagram ([3]) connects the ideal 
quantum computation, G, to the calculation of P, and, 
separately, we have that Diagram connects the ideal 
quantum computation, G, to the actual implementation 
of the quantum computation, P, that is realized by a 
practical physical device. To establish the utility of prac- 
tical quantum computers, what is needed now is a rela- 
tion that connects the actual, device-implemented quan- 
tum computation, P, to the intended calculation of F. 
By analogy with eq.((3]), we seek an expression of the form 



tr (Odm (P M (Idm (a:)))) = tr ^E F(x) P M (l dm (x)) ) > 1 - p' , (7) 

for some p' , where P M = Mj c ^i}PMji^ c }. This is established with the following calculation: 

tr(0 dm (P M (I dm (.x)))) = tr (y^^P^CW^))^^ ) 

= tr ^V F(x) {P M (l dm (x)) - Ul dm (x)U^}^E F(x) ) 
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+ ti (yjE F{x) UI din {x)U* yjE F{x) ) 
> -a + 1 — p, 

where we have used the property of the Schatten 1-norm, 



tr (yE F(x) {P M (I dm (z)) - Ul dm (x)tf}^E F(x) J < \\P M {l dm (x)) - C/I dm (a;)C/t||i , 
and the condition 

\\P M (l dm (x)) - C/Wz^lli = l|M {c ^ 1} (P • (M {1 ^ c} (I dm (.x)))) - G(I dm (z))||i < a 



(8) 



(9) 



(10) 



The physical meaning of this result can be transparently 
expressed by combining the content of Diagrams ^ and 
© to yield 



T(H, 
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This diagram is nearly commutative in the sense that 

tr(O dm (P M (I dm (a;)))) = tr(6 dm (P(I dm (x)))) 

> l-(p + a) , (12) 

where I d m = M{i^ c } o I dm , O dm = O dm o M{ c ^i}. We 
have thus established eq.0 with p' = p + a. Thus the 
total probability of failure that the overall computation 
yields the correct final result is bounded by the sum of 
two terms: (1) an upper bound on the failure probability 
for the overall computation when the quantum compu- 
tation is perfectly implemented without errors, and (2) 
an upper bound on the implementation inaccuracy of the 
actual quantum computation implemented by a practical 
device. This result is both intuitively simple and tech- 
nically subtle. It is not surprising that the total failure 
probability for the overall computation to yield the cor- 
rect final result is bounded by quantities that describe 
the two sources of error. However, it is not immediately 
obvious that these quantities should combine in such a 
simple way. On the one hand, the quantities p and p' rep- 
resent bounds on probabilities for obtaining certain out- 
comes from a measurement of the final quantum state. 
On the other hand, the quantity a is a bound, not on 
a probability, but on the normed difference between the 
states resulting from idealized and actual quantum com- 
putations, respectively. 

Diagram (Ti"Tj) and eq. (TT2"|) relate the actual implemen- 
tation P of the quantum computation to the intended cal- 
culation of F. They state the criterion for a realization, 
P, of a quantum computer, operating fault-tolerantly in 
the presence of residual errors, to correctly implement an 
instance of the overall computation. If p' is sufficiently 
small a majority vote algorithm will successfully identify 



the correct outcome F(x). For example, in the case Y is 
binary, the calculation of F succeeds by majority voting 
if < 1/2. 

Note that in the idealized limit in which error correc- 
tion perfectly and permanently removes all residual errors 
(i.e., in the limit a = 0), our result reduces to the corre- 
sponding result of the Kitaev model (i.e., eq.(TTJ]) reduces 
to eq.g])). 



IV. APPLICATIONS TO PRACTICAL 
SPECIFICATIONS FOR FAULT TOLERANCE 

A. General Result Relating Fault Tolerance 
Theory to the Overall Computation 

The above analysis, in which the constraints of fault 
tolerance are explicitly combined with those of the Ki- 
taev model, has important practical applications to the 
specification of error tolerances for quantum comput- 
ers. The designer of a quantum computer wishes to 
achieve some upper bound p' on the probability that 
the final measurement of the final quantum state pro- 
duced by the quantum computation fails to yield the 
correct final result for the overall computation. As dis- 
cussed above, there is some bound, p, on the inherent 
probability that the overall computation will fail even 
if the quantum computation is perfectly implemented. 
As discussed above, this is due to the intrinsic nondc- 
tcrminism inherent in the final measurement of the fi- 
nal quantum state resulting from the execution of the 
quantum computation, and it is an abstract property 
of the algorithm itself. Our result shows that the im- 
plementation of the algorithm by a realistic, i.e. noisy, 
quantum computer will meet the designer's success crite- 
rion provided that the implementation inaccuracy satis- 
fies||M{ c ^i } (P-(M { i_ +c} (p)))-G(p)||i < P '-p. Since the 
ideal failure probability bound, p, is a characteristic of 
the ideal quantum algorithm itself, this result effectively 
apportions the allowable noise in a quantum computation 
between a component, p, due to the inherent quantum 
mechanical indeterminacy associated to the measurement 
of the final quantum state, and a component, a, associ- 
ated to the dynamics of the quantum computation itself, 
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which includes other sources of noise, such as decoher- 
ence. Eq.© thus provides a success criterion for the 
design and implementation of a fault-tolerant quantum 
computation that, upon measurement of the final quan- 
tum state of which, produces the correct final result for 
the overall computation. 

As an example of how our result can be applied in or- 
der to achieve practical constraints on design parameters, 
suppose we wish to build a fault tolerant quantum com- 
puter which provides the correct solution to some specific 
problem (e.g. factoring a large number) with probability 
1 — p or better. That is, p is the required upper bound on 
the probability that the overall computation produces the 
wrong final result. From eq. (fT2|) we infer that the failure 
probability of the overall computation as implemented is 
bounded by p + a. We therefore require 



p + a = p 



(13) 



The quantity p is an intrinsic characteristic of the quan- 
tum algorithm, and can always, in principle, be deter- 
mined. We therefore require 



a =p-p , 



(14) 



which is a bound on the implementation inaccuracy of the 
quantum computation, sufficient to meet the prescribed 
success requirement for the overall computation. From 
eq.©, we see that this criterion is met if 

||M {C ^ 1} (P ■ (M {1 ^ c} (p))) - GOOId < P - P (15) 

for all p. 

Techniques of fault tolerance theory can be used to de- 
termine the probability £qc that the quantum computa- 
tion fails to produce the desired final quantum state. We 
therefore write the final state of the logical qubits result- 
ing from the actual, practically-implemented, quantum 
computation as 

M {c ^i } (P • (M {W (p))) = (1 - e QC )G(p) + e QC Perr , 

(16) 

where G(p) would be the result of an idealized quantum 
computation, p err arises from errors, and the entire right 
hand side of eq. (jTH]) is the state that results when errors 
occur. Eq. (p~5|) then becomes 



\\eQc\Perr-G(p)]\\l = e QC \\Perr ~ G(p) || 1 < p-p . (17) 

Since, 

\\perr ~ G(p)\\x <2 , (18) 

eqs. (TIT)) , (ri~5|) , and therefore (TT3"|) will be satisfied pro- 
vided 



1 ,„ 

£qc < ^ (P - P) 



(19) 



This is a new, quite general result relating the fault toler- 
ance theoretic parameter eqc to the constraints we have 



derived, which ensure that the overall computation yield 
the correct final result with some prescribed success prob- 
ability. This general inequality at once combines the con- 
straints dictated by the QCC, which apply to the dynam- 
ics of the quantum computation, with those coming from 
the Kitaev model, which apply to the measurement, in 
terms of the numerical parameters of standard fault tol- 
erance theory, which apply only to gate failure probabil- 
ities. 



Direct Example: Calculation of the 
Concatenation Level Function 



As an example of how our general result, Eq. p^)) . can 
be applied to the specification of practical fault-tolerance 
design parameters, we consider the fault tolerance ap- 
proach described in (4)], in which a concatenated quan- 
tum error correcting code is applied so as to reduce the 
failure probability at each level of concatenation. The 
failure probability for a single logical gate scales roughly 
as i 



£n — eth 



(20) 



where eo is the probability of failure for elementary gates, 
eth is the fault tolerance threshold, and is the failure 
probability of the gate at the Nth level of concatenation. 
For simplicity, in this example we assume that this re- 
lation is exact, that it applies equally to all gates, and 
that no other sources of error are present. If the quan- 
tum computation is comprised of N logical gates, then 
the failure probability for the quantum computation to 
yield the desired final quantum state is 



e QC 



Hi 



(21) 



Thus, in order to ensure that we obtain the correct final 
result to the overall computation with sufficient proba- 
bility of success, we make use of our general result in 
eq. tTT!)]) to require that 



< g CP - P) 



(22) 



is satisfied. From this we infer that the number of levels 
of concatenation sufficient to guarantee that the overall 
computation performs as required is given by: 



In 

N > log 2 p -=^- 



CO 



(23) 



This is a sufficient, not a necessary, condition. Alterna- 
tively, for a given level of concatenation, we could just as 
well derive a requirement on the error probability eo for 
elementary gates. In other words, this result enables us 



6 



15 



12.5 




2-10" 1 " 4-10" 10 6-10 10 8 - 10 10 1 - 10 



FIG. 2: Levels of concatenation, N, required to meet the 
performance criterion as a function of elementary gate error 
probability, eo. The performance criterion is to achieve the 
correct final result for the overall computation (not simply 
the desired final quantum state) with some prescribed success 
probability. In this example the success probability (p) is pre- 
scribed to be 0.6. The bound on error probability associated 
to the measurement following ideal quantum computation (p) 
is taken to be 0.2. The number of gates (K) is 10 12 , and the 
error threshold (eth) is taken to be 10 -9 . 



to explore the tradeoff between adding additional levels 
of concatenation and improving the performance of ele- 
mentary gates in order to achieve a given performance 
criterion in terms of the reliability of the final result. An 
example of such a tradeoff curve is shown in Figure [5] 
We stress that the points on the curve in Figure do not 
merely represent the amount of concatenation needed in 
order that the quantum computation produce a particu- 
lar final quantum state, but rather, fully incorporate the 
conditions that will ensure a successful overall computa- 
tion. 

We emphasize that the role of standard fault tolerance 
theory in our calculation is solely to express the error 
probability for the quantum computation in terms of el- 
ementary gate errors. Standard fault tolerance theory 
does not provide estimates of (nor bounds on) the im- 
plementation inaccuracy as denned in eq.fO, and, conse- 
quently, would not have allowed us to calculate the num- 
ber of levels of concatenation required in order to ensure 
that the overall computation produces the correct final 
result with the prescribed probability of success. 



V. CONCLUSION 

In this paper we have introduced a methodology for 
determining design parameters for a quantum computer, 
the values of which ensure that an overall computation 
of interest, comprised of an initial purely quantum com- 
putation followed by a measurement of the final quan- 
tum state produced by the quantum computation, yields 
a correct final result with some prescribed probability 
of success, as opposed to merely ensuring that the de- 



sired final quantum state is obtained. Thus, our method 
enables the determination of complete quantum compu- 
tational resource requirements associated to the actual 
solution of practical problems 

Our method fully accounts for two sources of error that 
may affect the final result of the overall computation: (1) 
errors which arise in the course of the quantum computa- 
tion itself, and (2) errors due to the indeterminacy intrin- 
sic to the final measurement which follows the quantum 
computation. Standard fault tolerance theory addresses 
the first type of error, but does not consider the sec- 
ond. Kitaev's model addresses the second type of error, 
but not the first. We have extended both standard fault 
tolerance theory and Kitaev's model, and have combined 
them, in order to yield quantitative results that reveal the 
interplay between the two. Although the analysis in this 
paper has been presented in the framework of the circuit 
paradigm for quantum computing, it is straightforward 
to apply our results to other paradigms, including the 
cluster state and adiabatic paradigms |8|. 

As a specific example of the practical application of 
our approach, we have explicitly calculated the number of 
levels of error correction concatenation needed to achieve 
a correct final result for the overall computation with 
some prescribed success probability. Extensions of the 
current calculation will include associating different fail- 
ure probabilities to different gates, as well as considering 
additional refinements dictated by imposing the QCC on 
the overall dynamics of the quantum computer [8| . 

Acknowledgments 

The authors thank F. Javier Thayer for his contribu- 
tions to the early phases of this research. They also thank 
Stephen P. Pappas and Anthony Donadio for their input. 
This research was supported under MITRE Technology 
Program Grant 07MSR205. 



7 



[1] P.W. Shor, SIAM Journal on Computing 26, 1484 (1997). 
[2] A. Kitaev, Russian Mathematical Surveys, 52, 1191 
(1997). 

[3] For simplicity we do not consider errors in the preparation 
of the initial quantum state in the present paper. 

[4] J. Preskill, Proc. Roy. Soc. Lond. A 454, 385 (1998). 

[5] D. A. Lidar, D. Bacon, K. B. Whaley, Phys. Rev. Lett. 
82, 4556 (1999). 

[6] M. Nielsen, I. Chuang, Quantum Computation and Quan- 
tum Information, Cambridge University Press, Cambridge 
(2000). 

[7] This can be illustrated with the example of the 7-qubit 
Steane code, for which H colnp is comprised of 7 compu- 
tational qubits for every logical qubit in illogical- We note 



that the linking maps M{i^ c } and M{ c _»i} do not represent 
transformations carried out by physical devices. The phys- 
ical encoding and decoding operations implied by the use 
of, e.g., the Steane code, must be performed by physical 
devices and are therefore subject to noise. Such physical 
encoding and decoding operations are properly included 
in P and not in the linking maps. In the context of the 
Steane code, one may think of the linking maps as rep- 
resenting the formal extension of Hilbert space to include 
the necessary ancilla qubits. 
[8] G. Gilbert, M. Hamrick, and Y.S. Weinstein, work in 
progress. 



